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Abstract 

We obtain static selfgravitating solitonic 3-brane solutions in the Einstein-Skyrme model in 
7D. These solitons correspond to a smooth version of the previously discussed cosmic p-brane 
solutions. We show how the energy momentum tensor of the Skyrme field is able to smooth out 
the singularities found in the thin wall approximation and falls fast enough with the distance from 
the core of the object so that asymptotically approaches the flat cosmic p-brane metric. 
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I. INTRODUCTION 



Models with extra dimensions have become intensively studied in the last few years, not 
only in string theory but in many other extensions of the standard model. Branes have 
played a major role on many of these models and have therefore been studied from several 
different perspectives Some of these branes have a description in terms of solitons in 

3j, monopoles 4], or even higher 



codimension objects jsl- One of the earliest models of particle physics to incorporate solitons 
to its spectrum, the Skyrme model 6|, has so far not been investigated in connection to the 
braneworld scenarios. The purpose of this letter is to study braneworld models based on 
the higher dimensional generalization of the Einstein-Skyrme model. 

One of the most interesting properties of the Skyrme model that makes it different from 
all the other topological defects mentioned above, is the presence of higher order corrections 
to the kinetic term. Recently there has been several studies on topological defects with 
non-canonical kinetic terms 7|]. Most of those models have been focusing on the differences 
introduced in the soliton solutions by the new term in the lagrangian. On the other hand, 
our model does not have a potential term in contrast to all the other cases studied in the 
literature of braneworld scenarios. Furthermore, in the absence of the Skyrme term, the 
theory does not present any smooth stable configuration even though there is a topological 
charge that one can define. It is only because of this new term that one is able to have a 
finite size stable configuration avoiding in this way the straighforward extension of Derrick's 
theorem 8|. It is then clear that the higher order terms in our Lagrangian are crucial for 
the solutions presented here and are not just a small correction to the model. 

On the other hand, the particular structure of the theory that we study here gives rise to 
solitonic solutions that are very much localized in the transverse directions. This property 
makes them good candicates to describe a smooth version of extended solutions previously 
studied only in the their thin wall limit. In fact, we will show that our thick brane solutions in 
the Einstein-Skyrme model asymptotically match the higher dimensional vacuum solutions 
of the pure Einstein's equations previously found in j^. 



the low energy theories as either domain walls {2I, strings 
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II. THE EINSTEIN-SKYRME MODEL IN 7D 



The action for our model is given by, 

Ses = [ d'X^g 



2h? 



R + C, 



(1) 



where R is the Ricci scalar, = 1/M|, with Mj denoting the 7-dimensional Planck mass, 
and Cs is the Skyrme Lagrangian density. 



,Tr{[LA,LB][L^,L'']), 



(2) 



where Fq and e are two free parameters of the model with units of [M]^/^ and [M] and 

La = U^OaU, (3) 

is the left chiral current and U G SU{2)} 

We are interested in finding the smooth solution for a 3-brane that is spherically sym- 
metric along the transverse directions in the bulk. We will also restrict ourselves to the four 
dimensional fiat brane solutions. Taking these constraints into account we can now write 
the most general metric of this form in the isotropic gauge as. 



where rj^y = diag{—l, 1, 1, 1). 



(4) 



We also impose the hedgehog ansatz for the chiral field, which is the natural spherically 
symmetric ansatz for Skyrme model; that is, we assume the following form for U{r), 



U{r) = cos/(r) + i i — ) r-' sin/(r) 



(5) 

where /(r) is the profile function to be solved for and with j = 1,2,3, are the Pauli 
matrices. 

Within this ansatz, the Lagrangian density then becomes, 

2 / n _:_,2 r \ 2 r / 2 



.F2 
-^0 



C^(r) \dr 



?i 11 + 



2 sin^ / 



sin^ / 



sin^ / 



62^2^2 (r)r2 



^ We use the following notation. The upper case latin indices A, B run over 0,..,6 and the greek indices 
II, v = 0, 3 denote the four dimensional spacetime coordinates. We use the the mostly positive signature 
and the Riemann tensor conventions of the form, Rbcd — '^c^bd ^ (^d^bc + ■■■■ 
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It is convenient at this point to rescale the radial coordinate eFor x and define, 



1 / , 2sin^/ 

u = „„ , , 1 



CHx) 



C^{x)x^ 



to obtain 



(6) 



Having simphfied the Lagrangian, the action for the Skyrme field becomes, 



Ss 



2nFo 



{uf + v)B\x)C%x)x'dxd^y. 



(7) 



It is now straightforward to obtain from this action the equation of motion for the field f{x) 
in the static case which is given by. 



where 



'B'{x) C'ix) u' 2 
4 — — + 3 — — H h - 

B{x) C{x) u X 



5uix) 



fix), 



(8) 



6v{x) 



(9) 
(10) 



and the primes denote the derivatives with respect to x. On the other hand, varying the 
Skyrme action with respect to the metric tensor qab yields 

2 SSs 



Tab 



where Tab is the energy-momentum tensor for the scalar field given by 

Tab = QabCs - ^Tt^LaLb) - ^ ^7^^Tr([LA, Lm] [L^, L^]). 



(11) 



In our ansatz it becomes, 
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V 

rpx 
X 

rriO 

rp<t> 
-L J, 



[uf - v) 



2C^{x) 

J-ei ■ 



-f 



12 



sin^ / 



(12) 



Using this energy- momentum tensor and our ansatz for the metric, (jl]), we obtain Einstein's 
equations of the form, 



3B"(x) 



2C"(x) ^ WHx) 



C'Hx) 



6B'(x) 



+ 



4C'{x) 3B'{x)C'{x) 



_B{x)C^{x) C^x) B^{x)C^{x) C^{x) B{x)C^{x)x C3{x)x 
{uf" + v) 5^^ 



8B'(x) 



QB'^(x) 



2C'{x) 8B'{x)C'{x) C'\x) 

B{x)C^{x)x ^ 52(x)C2(x) ^ xC^{x) ^ B{x)C3{x) ^ C^ix) 



^2 2p4 



Gn 



45" fx) 



+ 



G"(x) 



6B'Hx 



G"(x) 



AB'(x) 



G'{x) 



_B{x)G^{x) C3(x) 52(x)C2(x) G^{x) B{x)G^{x)x G'^{x)x 



K^e'F,' f , sin^ / 
2C2(x) V G^{x)x^ 



B{x)G'^{x) 



(13) 



Eqs. (fT3l) and ([8]) constitute the equations of motion for Einstein-Skyrme model consistent 
with the restrictions imposed by our ansatz. 

III. NUMERICAL RESULTS 



We want to find solutions for solitonic objects characterized by a topological charge Q 
that can be written in terms of an integral over the extra-dimensional space as. 



Q = ^ /" tr (L^L.L,) G\r)r'drdn^ 



TT 



/ Sin 2/ 
2 



/(oo) 



• (14) 



One can see that fixing the charge specifies the boundary conditions for our function /(r). 
In the following, we will consider the single charged soliton solution {Q = 1), which means 
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that we will take the following boundary conditions for the scalar field function /(r), 

/(0)=7r /(oo) = 0. (15) 



We want to integrate our equations of motion starting from the core of the defect, so we 
still need to specify the conditions for the metric functions at x = 0. We demand that our 
initial data at the origin does not have any singularity which in turn means that the most 
general expansion for the metric functions at x = should be of the form, 

B{x) = Bo + B2X^ + 0{x^) 
C{x) = Co + C2X^ + O(x^) 

f{x) = 'K + hx + hx^ + 0{x^). (16) 



Using this expansion in Einstein's equations, ( fT3l) . one can see that the higher order coeffi- 





cients can be obtained in terms of the parameters fi, Bq,Cq, k'^ = k^Fq and are given at the 



lowest order by, 

h 



/f(5Cg + ll/f)/l^ 
40Co 

ffil62ffk^ - 5Co^(-16 + 3k^) + 5Cg/f (8 + 9k^)) 



We numerically solve the system of equations ( 1131) and (E]) using the shooting method, 
i.e. adjusting /i, Bq and Cq such that the asymptotical solutions satisfy f{oo) = and 
B{oo) = C{oo) = 1. We show in Figs. 1-3, a sample of the numerical solutions found by 
this procedure.^ 



We found that similarly to what happens to the Einstein- Skyrme model in 4D [10|, 111 ] 
there is a critical value of k, k^, beyond which no more regular solutions exist. Our numerical 
investigation reveals that this critical value is around the order of k 



c2 ~ ^ 

c 20 



^ Note that we integrate the solutions to a much longer range in x where we clearly see the convergence 
of the different functions to their asymptotic values. We only plot a limited range in order to show the 
smooth structure of the soliton in the core region. 
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FIG. 1: Typical form of the B{x) function for the smooth 3-brane solutions. 
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FIG. 2: Typical form of the C{x) function for the smooth 3-brane solutions. 



Another interesting result is that for each value oi k < we have two branches of 
solutions. At k = kc, the two branches merge (See Fig. 4) such that beyond this point we 
always find a singularity at some value of x. This is another feature that is shared by the 



4:D system 



The numerical solutions we found asymtotically approach fiat space, and it is therefore 
possible to identify the form of the energy momentum tensor that sources these metrics, in 

n nn 

a similar way to what one does in AD spacetime [12[ . Following [13|, llJ] we obtain. 
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FIG. 3: Typical form of the f{x) function for the smooth 3-brane solutions. 
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FIG. 4: Two fundamental branches of soliton solutions. The shooting parameter /i is plotted as a 
function of k^. We see how the two solutions merge at k = kc- 



T^J''' =\im^(f nr],MK + ^^hi - d,h^) - d,h^J\ r^dn, (17) 

r '■oo J 

where Hab = Qab — Vab, denotes the deviation of our metric from fiat space, f* is the radial 
unit vector in the transverse 3-dimensional space and /i, i^, a = 0, ...,3 and i,j = 4,5,6. 
Using the expression of our ansatz we obtain in our case. 
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«2 




q-ADM 


fl 


Co 


Bo 




q-ADM 


flu 


Cqu 


Bqu 


0.01 


0.071614 


71.1456 


-2.21326 


1 


075 


0.9796 


0.120868 


120.077 


-35.9947 


4.007 


0.6524 


0.02 


0.139140 


69.1149 


-2.47453 


1 


165 


0.9550 


0.177366 


88.1029 


-21.2452 


3.506 


0.6802 


0.03 


0.202338 


67.0048 


-2.83941 


1 


283 


0.9299 


0.224255 


74.2627 


-14.3882 


3.070 


0.7091 


0.04 


0.260364 


64.6652 


-3.43789 


1 


460 


0.8940 


0.269732 


66.9919 


-9.92565 


2.633 


0.7450 


0.048 


0.303528 


62.8214 


-4.65254 


1 


766 


0.8445 


0.303671 


62.8510 


-6.52400 


2.140 


0.7952 


0.0485 


0.306311 


62.7438 


-4.88309 


1 


818 


0.8370 


0.306296 


62.7407 


-6.18143 


2.079 


0.8028 


0.0487 


0.307273 


62.6824 


-5.01605 


1 


847 


0.8329 


0.310717 


63.3849 


-6.02929 


2.055 


0.8060 


0.0488 


0.308737 


62.8519 


-5.10167 


1 


865 


0.8300 


0.308765 


62.8576 


-5.88885 


2.024 


0.8095 


0.0489 


0.308155 


62.6052 


-5.22771 


1 


892 


0.8268 


0.308099 


62.5938 


-5.74297 


1.996 


0.8131 



TABLE I: We show the values of the shooting parameters for different values of the gravitational 
coupling. We include both branches of solutions. The values of tq are given in C-Fq units and T^^^ 
in Fo/e units. 



rj^ADM ^ 8^ [r\W{r)B\r) + 2C(r)C"(r))] = -Tadm V,u (18) 

K, r— >oo 

It is then clear from this calculation that we can read off the value of the tension of 
the 3-brane source from the asymtotic behaviour of the metric. We show in Table I, our 
numerical results for the tension (Tadm) computed from the asymtotic form of the numerical 
functions B{x) and C(x), together with the values of the shooting parameters Bq, fi, Cq, 
for a range of k values. The subscript u denotes the upper branch in Fig. 4. 

The occurrence of the two branches in AD has been claimed to be linked to the existence 



of similar solutions in the Einstein- Yang-Mills system in 4D 15|]. It is likely that a similar 
situation may arise in our higher dimensional case. 



IV. COSMIC 3-BRANES IN THE ISOTROPIC GAUGE 

As we discussed in the introduction, our brane solutions have the same symmetry as 
the cosmic 3-brane gravity solutions discussed by Gregory 9|. Furthermore our branes are 
not charged with respect to any long range field, and therefore their energy momentum 
tensor is very well localized. It is then tempting to identify our solutions with the smooth 
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out version of the cosmic 3-brane examples. This possibihty was in fact aheady suggested 
by Gregory in j^. In the following we will prove that this is indeed the case, by showing 
that the metric solutions found numerically in the previous section match asymptotically 
the thin wall vacuum solutions found by Gregory. Let us briefly review the cosmic 3-brane 
geometries. The solutions of vacuum Einstein equations found in relevant for us are of 
the form, 



where 



ds^ = F{ff''r]^,dy>'dy'' + F{ff^df^ + f^F{ff^dnl (19) 



Fir) = 1 - ^ (20) 
r 



These analytic solutions are written in a different gauge from the numerical ones found 
in the previous section, but it is always possible to transform them into the isotropic gauge 
of the form, 

ds^ = B{rfr]^^dy^'dy'' + C{rY{dr^ + r'^dVt^) . (21) 
In this gauge, the 3-brane vacuum solutions become. 

Using the asymtotic form of this metric and the expression for the ADM energy-momentum 
tensor found in flTSl) we arrive at. 



T^r' = -^r,r^.. : (23) 

K 

which depends on the single parameter tq that completely characterizes the 3-brane metric 
solution. We give in Table 1 the results obtained for this parameter in our numerical 
examples. Note that our numerical solutions are smooth everywhere contrary to the analytic 
expressions in Eg. (1221) that have naked singularities located at r = tq and at r = 0. On the 
other hand, we will show that the analytic solutions are in fact a good approximation to our 
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FIG. 5: Comparison between the numerical B{x) and the thin wall solution (|22p . 
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FIG. 6: Comparison between the numerical C{x) and the thin wall solution ()22p . 

numerical results in the asymptotic region where r >> Tq but start deviating from them as 
r ~ tq. In order to see this, we proceed in the following way. 

We first identify, using Eq. (fT8|) and Eq. (123!) . the value of the parameter tq by looking at 
the asymtotic form of the ADM stress tensor of our numerical solution. Once this parameter 
is fixed, it singles out a particular member within the family of solutions given by Eq. (l22l) . 
Using this parameter we now plot in Figs. 5 and 6 the function B{x) and C{x) and compare 
them with the same functions in the numerical calculations. As we see the asymptotic form 
of the two functions, the numerical ones and the analytic vacuum solutions, are in perfect 
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agreement so we conclude that our numerical solutions do, in fact, match their thin wall 
counterparts. 



V. CONCLUSIONS 

We have constructed, using numerical techniques, static 3-brane solitonic solutions in the 
ID Skyrme model coupled to gravity. We have shown that their asymptotic form can be 
well approximated by the analytic vacuum solutions of pure Einstein theory obtained in . 
The presence of the energy momentum tensor that makes up the core of the brane is able 
to smooth out the singularities that show up in the analytic case, rendering these solutions 
completely smooth. 

The Skyrme branes obtained in this paper are asymptotically fiat and therefore represent 



a good candidate for regular braneworlds in the DGP model in ID |16l . Il7j ] . It is also clear 
that one could generalize our Skyrme model to higher dimensions in order to accomodate 
the DGP braneworld models of higher codimension within a similar framework to the one 
studied here. On the other hand, the properties of the gravitational sector of the braneworld 
may be affected by the details of the brane core [l4| so it would be interesting to test these 
ideas within our model. 

Furthermore, we have found that there is a maximum value of the coupling constant 
beyond which smooth solutions are not possible anymore. It is not clear what type of 
objects one should obtain for larger values of the coupling constant. One possibility is to 
relax the staticity of the metric. This seems to suggest that the branes would start to inflate 



in a similar way to the solutions found in 18|, ll9|, |20|. This is certainly a possibility although 



we note that the situation is slightly different from the usual defect solutions since we do 
not have any potential energy in our model, so the possibility of having topological inflation 
2l| at the core does not seem very likely. 

We have also found that there are two branches of solutions. We expect, as it happens in 
the 4D case, that the lower (upper) branch corresponds to a stable (unstable) configuration, 
although the stability calculation has to be performed taking into account the new possible 
channels opened due to the extradimensional nature of the solutions we have. 

We hope to come back to these and other interesting issues for Skyrme Branes in a future 
publication. 
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